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Five-component minimum-energy bound states and mobile vector solitons of a spin-orbit-coupled 
quasi-one-dimensional hyperline-spin-2 Bose-Einstein condensate are studied using the numerical 
solution and variational approximation of a mean-field model. Two distinct types of solutions with 
single-peak and multi-peak density distribution of the components are identified in different domains 
of interaction parameters. From an analysis of Galilean invariance and time-reversal symmetry of the 
Hamiltonian, we establish that vector solitons with multi-peak density distribution preserve time- 
reversal symmetry, but cannot propagate maintaining the shape of individual components. However, 
those with single-peak density distribution violate time-reversal symmetry of the Hamiltonian, but 
can propagate with a constant velocity maintaining the shape of individual components. 
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I. INTRODUCTION 

Bright solitons are localized wave packets that can tra¬ 
verse at a constant speed without changing their shape 
due to a balancing of nonlinear attraction and dispersive 
effects. Solitons have been studied in a variety of sys¬ 
tems which include water waves, non-linear optics, Bose- 
Einstein condensates (BECs), etc [T]. In case of BECs, 
the nonlinear attraction arises because of the net attrac¬ 
tive inter-atomic interactions, whereas the kinetic energy 
term in the Hamiltonian is the source of dispersion. In 
binary mixtures of scalar BECs, the presence of two com¬ 
ponents allows for the emergence bright vector solitons 
[2]. Solitons have also been studied in spinor BECs with 
hyperfine spin / = 1 [3] and / = 2 j4] in the absence of 
spin-orbit (SO) coupling. In a neutral atom, there is no 
coupling between the spin and the center-of-mass motion 
of the atom [S]. The SO coupling can be generated by 
controlling the atom-light interaction, and hence creating 
artificial Abelian and non-Abelian gauge potentials cou¬ 
pled to the atoms |3]. The first experimental realization 
of SO coupling [3, with equal Rashba [S] and Dressel- 
haus [3] strengths, in a neutral atomic BEC by dressing 
two atomic spin states with a pair of lasers paved the 
way for other experimental studies on SO-coupled spinor 
BECs jTU]. Recently, the effect of SO coupling on the 
solitonic structures has been studied in pseudospin-1/2 
UMa and spin-1 BECs la 

In this paper, we study an SO-coupled spin-2 BEC in 
a quasi-one-dimensional (quasi-lD) geometry character¬ 
ized by three interaction strengths cq oc (4 a2 + 804)/?, 
Cl (X (04 — 02)11, and C2 = (lao — 10a2 + 3a4)/7, where 
oo, 02, and 04 are s-wave scattering lengths in total spin 
= 0, 2, and 4 channels, respectively HSIIIZI- Our the¬ 
oretical analysis based on the SO-coupled single-particle 
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Hamiltonian allows us to identify two subdomains with 
distinct physical properties separated by the C 2 = 20ci 
line. In the C 2 < 20ci subdomain, the component den¬ 
sities of the wave function show a multi-peak structure, 
whereas in the C 2 > 20ci subdomain, a single-peak struc¬ 
ture of the same is found. Based on the solutions of the 
SO-coupled single-particle Hamiltonian, we identify an 
appropriate variational ansatz to calculate the minimum- 
energy bright soliton solutions for an SO-coupled BEC. 
The single-particle SO-coupled Hamiltonian is found to 
break Galilean invariance. By solving the Galilean trans¬ 
formed single-particle Hamiltonian dynamics, we demon¬ 
strate that only single-peak solitons can propagate with¬ 
out changing their shape. 

In order to find the moving vector solitons of the SO- 
coupled spinor BEC, we transform the SO-coupled GP 
equation using a Galilean transformation. Two degener¬ 
ate solutions of the GP equation in the rest frame evolve 
into two velocity-dependent nondegenerate solutions in 
the moving frame. We show that the multi-peak soliton 
of Ref. [13] is a linear combination of two such degen¬ 
erate solutions in the rest frame. In the moving frame 
the degeneracy is removed, and the multi-peak soliton 
cannot be formed by the same linear combination. The 
multi-peak soliton of Ref. M is thus velocity depen¬ 
dent and cannot propagate maintaining its shape. No 
such linear combination of two degenerate solutions in 
the rest frame is needed in the formation of the single¬ 
peak soliton identified in this paper. These solitons can 
emerge as the ground states of the coupled GP equation 
in a moving frame, and hence can move with a constant 
velocity maintaining their shape, including the densities 
and phase profiles of individual components, unlike the 
moving multi-peak solitons of Refs. [121114] which show 
spin-mixing dynamics. 

In Sec. 1^ we describe the mean-field formalism used 
to study the SO-coupled spin-2 BECs. Here, taking 
into account the effect of interactions on the solutions 
of the non-interacting and trapless system, we identify 
two types of ground state solutions with single-peak and 
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multi-peak densities, respectively. In case of ground-state 
solutions with single-peak densities, the reduction of the 
conpled GP eqnation to a single nonlinear partial differ¬ 
ential equation nsing single-mode approximation (SMA) 
is also discnssed in this section. In SMA all the compo¬ 
nent wavefnnctions are assnmed to have the same spatial 
dependence but can have different normalization |15| . In 
Sec. |III[ we study the minimum-energy stationary soli- 
tons by variationally minimizing the energy. We also 
investigate the moving solitons from a consideration of 
Galilean invariance of the Hamiltonian. In Sec. EYi we 
present the nnmerical resnlts for both types of solitons. 
We conclude by providing a snmmary of the stndy in Sec. 

El 


II. SPIN-ORBIT COUPLED BEC IN A 
QUASI-ID TRAP 


Although, unlike in the case of charged electrons in 
an atom, an SO conpling interaction does not naturally 
arise in the case of nentral atoms in a spinor BEG, an SO 
conpling can be realized by an engineering of atom-light 
interaction. An SO-coupled spinor BEG in a quasi-lD 
trap can be realized by making the trapping freqnencies 
along the two axes, say y and z (wyjWz), much larger 
than that along the x axis (wa,). This leads to a BEG 
in which the dynamics is restricted only along x axis. 
The single-particle Hamiltonian of this BEG with in a 
qnasi-lD trap is mm 

Ho = ^+ V{x) + (1) 


where Px = —ihd/dx is the momentnm operator along x 
axis, V(x) = mujxx'^ f2 is the harmonic trapping potential 
along x axis, and is the irreducible representation of 
X component of the spin operator: 
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X 


/O 1 0 0 0 \ 

1 ^ 0 0 

0 0 a/372 0 

0 0 0 1 

V 0 0 0 1 0 / 


( 2 ) 


level and the conpling was achieved after introducing the 
nonlinear spinor interactions in the BEG. 

Taking interactions in the Hartree approximation and 
nsing the single-particle Hamiltonian Q, a quasi-ID 
spin-2 BEG can be described by the following set of five 
conpled mean-field partial differential eqnations for the 
wave-function components m 
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+ j ++ c+^'i} + (5) 


where cq = 2fi^(4o2 -|- ci = 2h?{ai — 

a2)/{lmll^), C2 = 2h^{7ao - 10a2 -f 3a4)/{7mllJ), Oq, 
02 and 04 are the s-wave scattering lengths in the to¬ 
tal spin /tot =0,2 and 4 channels, respectively, Pj{x) = 
\ipj\‘^ with j = 2,1,0,-1,-2 are the component densi¬ 
ties, p{x) = Y^j =-2 Pi total density, and = 

•^Jh/(mujyz) with ujyz = y/ujyOJz is the oscillator length in 
the transverse y — z plane and 

F+=F*_=FxF iFy 

= 2('02'0i + ip-iip-2) + V 6 {iptipo + ipoip-i), 

Fz=2{\lP2?-\lP.2?) + \lP,\^-\lP.,\\ 

_ 2lp2lp-2 - ‘ 2 . 1 pl 1 p-l + ipo 
V 5 


The SO conpling corresponds to eqnal strengths of 
Rashba {HxPx + '^yPy) Dresselhaus (F,xPx ~ "^yPy) 
couplings [TBl - ED] . where Hy is the irreducible representa¬ 
tion of y component of the spin operator. An equivalence 
of these couplings with the usnal Rashba (T,xPy — YiyPx) 
[S] and Dresselhans {TixPy -I- '^yPx) [S] conplings can be 
shown by a rotation in the spin space HHiiini. 

In onr previons studies ED El, we employed a dif¬ 
ferent SO coupling (px'^z) again with equal Rashba [S] 
and Dresselhans d strengths. The present single-particle 
Hamiltonian Q already conples the different spin compo¬ 
nents. The previons SO interaction EDUa did not cou¬ 
ple the different spin components at the single-particle 


where Fx,Fy, Fz are the three components of the spin- 
density vector F, and 0 is the spin-singlet pair amplitude 
|20| . Before proceeding further, let us transform the Eqs. 
into dimensionless form using 

i=uJxt, i= (Pj{x,i) = ^^ipj{x,i), (6) 

to V-/V 

where Iq = ^/h/JjnuPx) is the oscillator length along x 
axis, and N is the total number of atoms. Then, Eqs. 
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in dimensionless form become 


.d(j)±2 ( I A ^ 

+ Ci(i?^(/)±l ±2F^(/.±2) + {C2/V5)er^2, 
. 50±1 f 1 92 .V 


d(t>±i 

dx 

d(j>±2 

dx 


(7) 


- * 7 "^^ + Cl (vW2F^</>o + F±cji±2 ± 

- (52/v^)0<^;i, ( 8 ) 



.i9(?!)o 
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\ 2 9x2 


+ 17(i) + Cop 


</>o 




H+ ^ci(F_(^_i + F^0i) + ^0(/)g,, (9) 


where F = i^/ 2 , 7 = hkr/[muixlo), cq = 2 iV( 4 a 2 + 
3a4)lo/{7llJ, Cl = 2iV(a4 - a2)hlill'll), C 2 = 
2iV(7ao - IO 02 + 3a4)lol{Jll^), Pj{x) = \(j)j\‘^ with j = 
2 , 1 , 0 ,- 1 ,- 2 , and p{i) = Ei =-2 and 

= FI =2ir2^4 + Cl<^- 2 ) + 

( 10 ) 

F,-=2(|<(>2p-|<^-2p) + |0 ip-|0 -iP, (11) 

0 _2 ^2<(’-2 ~ 2(()i(()-i + (^Q 

The total density is now normalized to unity, i.e., 
p{x)dx = 1. We will represent the dimensionless 
variables without tildes in the rest of the manuscript for 
notational simplicity. 

For a non-interacting system in the absence of a trap¬ 
ping potential, there are five linearly independent solu¬ 
tions of Eqs. 0 -@: 


$1 = 

ikx T 

— ( 1 , 2 , 76 , 2 , 1 ] , 

(13) 

$2 = 

ikx T 

— ( 1 ,- 2 , 76 ,- 2 , 1 ] , 

(14) 

$3 = 

pikx 

— (- 1 ,- 1 , 0 , 1 , 1 )^, 

(15) 

$4 = 

Akx 

— (- 1 , 1 , 0 ,-!, 1 )^, 

(16) 

$5 = 


(17) 


where T stands for transpose. Here $1 and $2 are two 
degenerate solutions corresponding to energy per particle 
E/N = k‘^/2 ± 2kj, which has a minimum E-aan/N = 
— 27 ^ at fc = t 27 , where the upper and lower signs are 
for $1 and <i> 2 , respectively. Similarly, the eigen functions 


FIG. 1: (Color online) Energy dispersion E/{N'y^) vs. fc/y for 
the eigen functions of the single-particle SO-coupled Hamil¬ 
tonian for the five wave functions >l>i, where i = 1,..., 5. 


$3 and $4 are also degenerate with energy E/N = fc^/2± 
fcy, which has a minimum E^-m/N = — 7^/2 at fc = =P 7 
with upper and lower signs for $3 and $ 4 , respectively; 
whereas the non-degenerate eigen function $5 has energy 
E/N = fc^/2, which has a minimum E-^ain/N = 0 at 
fc = 0. The E/{Nj'^) vs. fc /7 dispersion curves are 
shown in Fig. ^ and we will consider the lowest-energy 
ground state $12 with fc = = 727 . 

Compared to the SO-coupled pseudospin-1/2 p5l [21] 
and spin -1 condensates [23], there are two sets of degen¬ 
erate eigen functions in the spin-2 case. Due to this, two 
distinct types of muti-peak stationary profiles can emerge 
in SO coupled spin-2 condensate. 

The most general solution of Eqs. 0 -® with mini¬ 
mum energy for a BEC with a uniform density n in the 
absence of interactions and trapping potential can be ob¬ 
tained by the linear superposition of ^-nd yEi ^2 
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(19) 




where |aip -I- |a 2 p = 1 , so that $ is normalized to 


unity. The solution (191 implies that the magnetization 


A4 = J F^dz = 0. It implies that we can only obtain sta¬ 
ble bright solitons with zero magnetization for the SO- 
coupled spinor condensate considered in this paper. On 
the other hand, in the absence of SO coupling, one can 
have stable bright solitons with non-zero magnetization 
for spin -2 condensate jj). 

The time-reversal-symmetric single-particle Hamilto¬ 
nian 0 violates parity. This will have interesting conse¬ 
quences on its eigen functions. The effect of time-reversal 
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symmetry operator T acting on a spin or orbital angular- 
momentum state Ij, m) is |2()| 

( 20 ) 


Time-reversal symmetry of a quantum state 4* requires 
that it should be the same as its time-reversed state 
apart from a phase factor. The degenerate states (13) 
and (141 violate time-reversal symmetry and are con¬ 
nected by time-reversal symmetry operator: T|$i) = 
— |$ 2 ),T|$ 2 ) = —l^hi)- In general, for arbitrary ai and 
02 the states (19) do not have time-reversal symmetry. 

In the presence of interactions, the interaction energy 
per particle is j^D] 




Co , „ 
= —n + 2n 


Cl 


|2 I |2 ( C2 20ci \ 

|ai| |a 2 | 1 -^- 1 


.( 21 ) 


If C 2 < 20ci, the minimum of ejnt corresponds to \ai\ = 
|a 2 | = l/v^- This state is nondegenerate and has multi¬ 
peak density distribution for the wave-function compo¬ 
nents and is time-reversal invariant. The spin expec¬ 
tation and absolute value of spin-singlet pair amplitude 
per particle for this state are, respectively, |F|/n = 0 
and |0|/n = l/'Jb. On the other hand, for C 2 > 20ci, 
the Cint can be minimized if |ai| = 1 , |q; 2 | = 0 or 
I Oil = 0,1 02 1 = 1. These two plane-wave states are de¬ 
generate and are connected by the time-reversal operator 
and hence break time-reversal symmetry of the Hamil¬ 
tonian. The spin expectation per particle for this state 
|F|/n = 2, and the absolute value of spin-singlet pair am¬ 
plitude per particle |0|/n = 0. These states have single¬ 
peak density distribution for the wave-function compo¬ 
nents and can also be studied using SMA. In SMA, the 
order parameter can be written as 


= -(1,2, Ve, 2, l)'^())sMA(a:,t), (22) 


based on the minimum-energy solutions (131 of the 


single-particle Hamiltonian. For SMA to be valid, all 
the spin-dependent interactions should be much smaller 
than the spin-independent interactions m- Hence, us¬ 
ing Eq. (22) in Eqs. (|^-(|^ and neglecting the ci- and 
C 2 -dependent terms, one can obtain the single nonlinear 
differential equation 


. i9<^sma 
* dt 


1 >9^ T,. N 5 ' 


<)'sma,(23) 


where the tildes have been dropped, 
be solved numerically. 


Equation (23) will 


III. BRIGHT SOLITONS 
A. Variational analysis 

We just found that the physical properties of the sys¬ 
tem are different for C 2 > 20ci and C 2 < 20ci, and now we 


present a variational analysis of the bright soliton with 
minimum energy in these two domains. We variationally 
minimize the energy of the trapless BEG Ennsg, given 
by 


E = N 



■ . d4>-l 

. i*d(j)2 . d(j)-2 

-.70,—-27^-,^ 


(24) 


in the two domains which are separated by the C 2 = 20ci 
line. An appropriate variational ansatz 4>var is con¬ 
structed by taking the product of the superposition of 
the eigen functions corresponding to minimum energy in 
Eqs. (13) and (14) with a localized spatial soliton, i.e.. 


$var = — 4’sech(CTx), 


(25) 


where cr is the variational parameter characterizing the 
width and strength of the soliton, and |aip -|- |Q; 2 p = 1- 
As discussed in Sec. |TTj the exact values of ai and a 2 
differ in the two domains. In order to search for the stable 
stationary bright solitons with higher energies, one can 
also consider the following variational ansatz: 


< 1 > = 
^ va.r 


^ var 

l^var 


-(q!i$3 -I- a24)4)sech(crx). 


4>var = -^<i)5sech((Tx). 


(26) 

(27) 


Ansatz (26) and (27) can lead to bright solitons with suc¬ 


cessively higher energies than the one described by Eq. 
(25). Nevertheless, in the present paper, our emphasis is 


to calculate the minimum-energy bright solitons consis¬ 
tent with Eq. (25). 

If C 2 < 20ci, we choose |ai| = |q! 2 | = l/-\/2 in Eq. 


(25). Then, the energy of the soliton is 


E = — (— 667 ^ -I- 5cr^ -I- 5crco -I- (TC 2 ) 

oU 


(28) 


with a minima at 


CT = (5co -kC 2 ) 


(29) 


provided that 5co -I- C 2 < 0. The energy and hence the 
shape of the soliton are not the functions of interaction 
parameter ci due the fact F = 0 for |ai| = \a 2 \ = l\/2 as 
has been discussed in Sec. [Hj This parameter domain cor¬ 
responds to the antiferromagnetic phase (F = 0,]^/n = 
lA^) El. The state $var, defined by Eqs. ( |25| and 


(19) with |ai| = |q( 2 | = l/'sA is time-reversal symmetric. 
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If C 2 > 20ci, we choose |ai| = 1, |a 2 | = 0 or vice versa 
in Eq. (25). The energy of the soliton, then, is 

E = ^ (- 127 ^ + cr^ + crco + 4crci) (30) 
with a minima at 


CT = -- (co + 4ci), 


(31) 


provided that cq + 4ci < 0. The shape of the soliton 
is independent of interaction parameter C 2 due vanishing 
0. This parameter domain corresponds to the ferromag¬ 
netic phase (|F|/n = 2, 0 = 0) [T7]. The $var with 
|ai| = 1 , |a 2 | =0 or vice versa leads to single-peak 
density distribution for wave-function components and 
breaks time-reversal symmetry. 


B. Moving bright solitons 


The GP equation governing the dynamics and statics 
of a scalar BEG is Galilean invariant. The implication 
of this invariance can be understood by considering the 
scalar ID GP equation in dimensionless form 


.d(j)(x, t) 
* dt 


1 d'^(j){x, t) 

2 dx'^ 


c\(l>{x,t)\'^4>{x,t), 


(32) 


where the nonlinearity is considered attractive (c > 0 ) 
for the formation of a bright soliton. This equation has 
the analytic solution 


cj){x,t) = x/(^sech [{x - vt)^] ( 33 ) 


where v is the velocity and a represents the width and 
the strength of the soliton. It implies that the stationary 
soliton (j){x, 0) = (cr/c)sech [x-ya] moves as the bright 
soliton (j)M{x,t) defined by [T^] 


^Mix, t) = (l){x - vt, 


(34) 


maintaining the width and strength (cr) fixed. In the 
case of a multi-component spinor BEG, Eq. (34) remains 


valid with a multi-component wave function $ replacing 
the single-component (j) while the multiplying exponential 
factor remains unchanged. 

Now, let us examine the Galilean invariance of the SO- 


coupled Hamiltonian. Equation (34) implies that for the 


Galilean transformation x' = x + vt, t' = t, where v is 
the relative velocity of unprimed coordinate system with 
respect to primed coordinate system, the wave function 
4) of Eqs. (|^-(|^ should transform to 4 >m related by 

4>(x, t) = $m(x , ^)*72_ ( 35 ) 


Substituting Eq. (35) in Eqs. 0 -@ and using d/dx = 
djdx', djdt = 9/i9F + vdjdx', we get 


.d^uix' ,t') 
* dt' 


■ 1 92 
2 9a:'2 





^Mix',t'), 

(36) 


where, for the sake of simplicity, we have suppressed the 
terms proportional to cq, ci, and C2, which should remain 
unchanged. We have also dropped a a dependent addi¬ 
tive term in the Hamiltonian which does not contribute 
to the dynamics. 

The extra term on the right hand side of 

the equation indicates that the SO-coupled Hamiltonian 
is no longer Galilean invariant. The SO-coupled equation 
( [M| has the solutions $1 and 4)2 of Eqs. (13) and (14) 
with energies E = —2N{v^ + 7 ^) and E = 2N{v^ — 
respectively. Hence, the degeneracy between 4)i and 4>2 
is removed for v ^ 0. This in turn implies that the 
superposition of 4?! and ^2 in Eq. (18) in the rest frame, 
which leads to a multi-peak soliton, is not possible in the 
moving frame. In other words, the multi-peak soliton 
cannot propagate with a constant velocity maintaining 
its shape. On the other hand, for a single-peak soliton no 
mixing between 4>i and $2 is allowed based on energetic 
considerations as discussed in Sec. Ulland Sec. IHII Hence 
a single-peak soliton can traverse with a constant velocity 
maintaining its shape. 


C. Numerical solutions 


We numerically solve the coupled Eqs. ([^-([^ using 
split-time-step Grank-Nicolson method |26l [27] in imagi¬ 
nary and real times. The ground state is determined by 
solving Eqs. 0-0 in imaginary time. In order to fix 
both the magnetization and normalization, we use the 
approach discussed in Ref. m- Accordingly, after each 
iteration in imaginary time, we renormalize the compo¬ 
nent wave functions as 


(j)j{x,T + St) = dj(j)j{x,T), 

(37) 

where dj’s satisfy the following relations 

did-i — dl, 

(38) 

d2d-2 = do, 

(39) 

d 2 (i^i = dl, 

(40) 

and 

d\N2+dldlNi + dld\No + dldlN.i + dgiV _2 

= 1, (41) 

2dlN2 + dldlNi - dldlN_i - 2d^V_2 = M. 

(42) 


Hence, one needs to solve the coupled set of non-linear 
Eqs. (41 )-(|42|) to determine do and di, which can be back 


substituted in Eqs. (38)-(40) to calculate the remaining 
normalization constants. The spatial and time steps em¬ 
ployed to generate the numerical results in this paper are 
Sx = 0.05 and 6t = 0.0000625, respectively. 


IV. NUMERICAL RESULTS 

First, we consider an SO-coupled / = 2 spinor BEG 
of 10000 23Na atoms trapped in a harmonic trapping po¬ 
tential with u}x/{2,n) = 20 Hz, ojyl{2Tr) = uJzI{2'k) = 400 
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FIG. 2: (Color online) Ground state density profiles of an 
SO-coupled trapped spinor BEG for (a) co = 242.97, ci = 
12.06, C 2 = —13.03 < 20ci, and (b) co = 201.36, ci = 
-1.8,C 2 = 24.15 > 20ci. 


Hz. The oscillator lengths with these parameters are 
Zo = 4.69 and ly^ = 1.05 /im. This value of Zg has been 
used in all calculations to write the dimensionless GP 
equations 0 -§. The scattering lengths of ^^Na in total 
spin /tot = 0,2 and 4 channels are og = 34.90^,02 = 
45 . 803,04 = 64.503 im HO], respectively, resulting in 
cg = 242.97, Cl = 12.06, C 2 = —13.03 < 20ci. Here 03 is 
the Bohr radius. The ground state density profile with 
A4 = 0 is shown in Fig. [^a). The multi-peak nature 
of the solution, obtained as result of the superposition of 
two counter-propagating plane waves, is consistent with 
analytic results obtained in Sec. |TT] The solution is dy¬ 
namically stable and retains its density and phase pro¬ 
file if real-time propagation of the dynamics is performed 
upon small perturbation. 

Next we consider the trapped BEG in the C 2 > 20ci do¬ 
main. For this, we consider og = 52 . 3503,02 = 45.8, and 
04 = 43 o 3 , which results in cg = 201.36, ci = —1.8, C 2 = 
24.15 > 20ci. The necessary modification of the scatter¬ 
ing lengths can be achieved by the optical and magnetic 
Feshbach resonance techniques |2S]. The trapping po¬ 
tential parameters and number of atoms are the same 
as those in Fig. [^a). The ground state densities with 
A4 = 0 are shown in Fig. I^b). The component den¬ 
sities obtained by using SMA) i.e., by solving Eq. (231, 
are in good agreement with the full numerical solution of 



X 
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FIG. 3: (Golor online) Numerical (num.) and variational 
(var.) densities of an SO-coupled / = 2 BEG soliton for 
(a) Co = —1.55, Cl = 0.39, C 2 = 0.48, (b) co = 4.94, ci = 
— 1.80,C2 = —8.06, and (c) co = — 2.74, Ci = —0.16, and 
C 2 = 1.77 > 20ci. In (c) the densities obtained by SMA 
are also shown. 


Eqs. Q-([^. The solutions in this case are dynamically 
stable plane wave solutions, resulting in a single-peak 
density distribution for the wave-function components, 
consistent with the analytical analysis in Sec. |TT] 

In order to obtain the bright solitons {V{x) = 0) in 
SO-coupled spinor BEGs, we consider the two cases high¬ 
lighted in Sec. HI A (a) C 2 < 20ci and Scg -I-C 2 < 0, and 
(b) C 2 > 20ci and Cg -I- 4ci < 0. In case (a), we consider 
og = —1.503, 02 = — 1 . 203 , 04 = 0, and N = 5000, 
which results in cg = —1.55,ci = 0.39, C 2 = 0.48. The 
numerical results for the component densities of the soli¬ 
ton are illustrated in Fig. [^a) together with the corre¬ 
sponding variational results, defined by Eqs. (251 and 
(291, with M. = Q and |ai| = \a 2 \ = l/'x/S- The solution 
in this case has multi-peak density distribution, is time- 
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FIG. 4: (Color online) (a) Nonzero numerical densities of an 
SO-coupled / = 2 BEG soliton with energy larger than the 
bright soliton shown in Fig. §a) . (b) Nonzero numerical 
densities for an / = 2 BEG soliton in cyclic phase in the 
absence of SO coupling. The interaction parameters for both 
(a) and (b) are cq = —1.55,ci = 0.39, C2 = 0.48. 


reversal symmetric, and dynamically stable. Despite the 
modulation in densities of the individual components, the 
total density profile (p) of the bright soliton as shown 
in Fig. j^a) does not have any modulation unlike the 
bright solitons discussed in Ref. HU. In this sense the 
present multi-peak solitons are spin-2 analogues of the 
stripe phase discussed in Ref. |23[ Moreover, in 

Ref. im the authors noted component symmetries like 
Rea^lk^) = — Real('I'^), and Imag(lE'^) = Imag('I'^), for 
the real and imaginary parts, which is not possible in the 
present model. This is because of different SO coupling 
used in the two studies, e.g., 'ypx'^x in this paper and 
'yPxCTz in Ref. HU. Had we used the SO coupling jpx'Sz 
the component symmetries of Ref. m can be obtained 
as noted in Table 1 of Ref. HU • 

In order to obtain the stationary bright soliton with 
higher energy, we use imaginary time propagation with 
($3 -p $4) exp(—a:^/2)/(-\/2-y7r) as the initial guess for 
the order parameter. The dynamically stable bright soli¬ 
ton thus obtained with the parameters of Fig. [^a) 
is shown in Fig. |^a). The period of density modu¬ 
lation in this case is twice of the period in Fig. |^a). 
Hence, there are two distinct stable multi-peak solitons 
in SO-coupled spin-2 condensate as compared to SO- 
coupled spin-1 condensate which can have only one type 
of multi-peak soliton HU. There also exists a stable 
single-peak soliton, not shown here, with order parameter 
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FIG. 5: (Golor online) (a) Dynamics of a soliton of Fig. 
I^a) placed at a; = —20 and set into motion with a con¬ 
stant velocity by multiplying the wave function by a phase 
Here co = -1.55, ci = 0.39, C 2 = 0.48, M = 0, 
and V = 1. (b) Elastic collision between the two moving 

solitons of Fig. ib) placed at a; = ±20 and set into mo¬ 
tion by multiplying by phase factors respectively. Here 

Co = 4.94, Cl = —1.80,C2 = —8.06, JV = 5000, A4 = 0, and 
V — 1. 


^/3p(x)/8(l, 0,—^/2/3, 0,1) and energy higher than the 
aforementioned two multi-peak solitons. For the same set 
of parameters, these three types of solitons — minimum- 
energy ground and higher-energy excited states — have 
the same total density p(x). 


In the absence of SO coupling the minimum-energy 
soliton corresponds to time-reversal symmetry breaking 
cyclic phase (|F| = |0| = 0) for ci > 0 and C2 > 0, 
C2 < 20ci |17L 1201121] . This is shown in Fig. Sb) for 
the interaction parameters corresponding to Fig. ia) 
and 7 = 0. The presence of SO coupling results in 
the time-reversal symmetric antiferromagnetic phase for 
which |F| = 0 and |0| > 0. Next we consider an 
SO-coupled soliton with oq = 3.das, 02 = 4.6805, 
04 = —l.Oas, and JV = 5000, which is equivalent to 
Co = 4.94, Cl = —1.80,C2 = —8.06, C2 > 20ci. The nu¬ 
merical results in this case are contrasted in Fig. .ib) 
with the variational results, defined by Eqs. ( [2^ ) and 
(311 with |ai| = 1, |a2| = 0 or vice versa. The soli- 
tonic solution here breaks the time-reversal symmetry of 
the Hamiltonian, as can be seen, for example, from Eqs. 
(251 and (131 with 02 = 0. In the absence of SO coupling 
the density distribution remains unchanged in this case 
for the bright soliton with zero magnetization, but the 




































component wavefunctions have constant phase instead of 
constant phase gradient. The solution thus retains the 
time-reversal symmetry of the Hamiltonian. The SMA 
can not be applied in this case to calculate the bright 
solitonic solution as cq > 0. To study the applicabil¬ 
ity of SMA, we consider cq = —2.74, ci = —0.16, and 
C2 = 1.77 > 20ci. The bright soliton obtained by the nu¬ 
merical solution of Eqs. ([^- (joL their variational approx¬ 
imation, viz. Eqs. ( |25[) and ( |31[ ) with |ai| = 1, \a 2 \ = 0, 
and SMA, viz. Eq. ( |23[ are shown in Fig. Sc). 

In order to set the soliton into motion with a constant 
velocity v, we multiply the static solution by exp {ivx) 
with V = 1. We find that bright solitons with standing 
wave density profile in the C2 < 20ci and 5co -I- C2 < 0 
parameter domain, exhibit spin-mixing dynamics as is 
shown in Fig. ^a) for the soliton with same interaction 
parameters as in Fig. ^a) and moving with velocity v = 
1 in dimensionless units. Thus, this type of soliton does 
not strictly preserve its shape while in motion. On the 
other hand, moving solitons obtained by multiplying the 
plane wave solutions in the cq + 4ci < 0 and C 2 > 20 ci 
domain, by exp(±iua;) strictly preserve their shape. This 
is evident from Fig. ib) , where we numerically study the 
collision between two bright solitons, initially located at 
X = ±20, and moving in opposite directions with a speed 
u = 1 (in dimensionless units). The solitons collide at 
a; = 0 and pass through each other without suffering a 
change in their shapes as is evident from Fig. ib). The 
interaction parameters in this case are the same as in Fig. 
ib). As discussed in the Sec. |IIIB[ these moving five- 
component vector solitons do not exhibit any spin-mixing 


dynamics. 


V. SUMMARY 


We study the generation and propagation of five- 
component vector solitons in an SO-coupled spinor BEG 
with hyperfine spin / = 2 using a mean-field GP equa¬ 
tion with three interaction strengths: cq, ci, and C2. Two 
types of solitons — single-peak and multi-peak — emerge 
in this case for C2 > 20ci and C2 < 20ci, respectively. In 
the former case, the solutions violate time-reversal sym¬ 
metry, whereas in the latter case, the solutions are time- 
reversal symmetric. The GP equation for this system is 
demonstrated to violate Galelian invariance. Analyzing 
the Galelian invariance of this equation, we show that 
the single-peak solitons can move with a constant veloc¬ 
ity maintaining constant component densities. On the 
other hand, the multi-peak SO-coupled solitons change 
the component densities during motion. 
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